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ABSTRACT 

Recent exam of large samples of ui Cen giants shows that it shares with mono-metallic globular 
clusters the presence of the sodium versus oxygen anticorrelation, within each subset of stars with iron 
content in the range -1.9 ^ [Fe/H] ^ -1.3. These findings suggest that, while the second generation 
formation history in oj Cen is more complex than that of mono-metallic clusters, it shares some key 
steps with those simpler cluster. In addition, the giants in the range -1.3<[Fe/II] -0.7 show a direct 
0-Na correlation, at moderately low O, but Na up to 20 times solar. These peculiar Na abundances 
are not shared by stars in other environments often assumed to undergo a similar chemical evolution, 
such as in the field of the Sagittarius dwarf galaxy. These O and Na abundances match well the yields 
of the massive asymptotic giant branch stars (AGB) in the same range of metallicity, suggesting that 
the stars at [Fe/H]>-1.3 in u Cen are likely to have formed directly from the pure ejecta of massive 
AGBs of the same metallicities. This is possible if the massive AGBs of [Fe/H]>-1.3 in the progenitor 
system evolve when all the pristine gas surrounding the cluster has been exhausted by the previous 
star formation events, or the proto-cluster interaction with the Galaxy caused the loss of a significant 
fraction of its mass, or of its dark matter halo, and the supernova ejecta have been able to clear the 
gas out of the system. The absence of dilution in the metal richer populations lends further support to 
a scenario of the formation of second generation stars in cooling flows from massive AGB progenitors. 
We suggest that the entire formation of lu Cen took place in a few lO^yr, and discuss the problem of 
a prompt formation of s-process elements. 

Subject headings: globular clusters: general — globular clusters: individual(a; Cen) 



1. INTRODUCTION 

A number of surprising observational findings of 
the latest dozen years suggest a completely new 
paradigm for the formation of Globular Cluster (GC) 
stars. The abundances of stars in most Galactic 
GCs observed so far show particular patterns, such 
as the oxygen vs. sodium and magnesium vs. alu- 
minum anticorrelations, not shown by the halo stars 
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(jCratton et al.l [20041), both in giants and in scarcely 
evolved stars (e .g. Grattoi Tet al.ll200lHRamirez fc Cohenl 
[20M iGTattoTet al.., 20041 : iCarretta et al.ll2009al|b[ ). The 
observed amount of chemically anomalous stars and the 
relative number of normal and anomalous stars suggest 
a sequence of events starting with the birth of a first 
stellar generation (FG) in a "protocluster" much more 
massive than todays' cluster, followed by the birth of a 
second generation (SG) from gas having the same iron 
content, but including hot-CNO cycled matter ejected 
by the FG stars evolving during the firs t phases of the 
clust er life (for a summary, see, e.g. iD'Ercole et al.l 
|2008[ ). Theoretical studies have identified the matter 
constituting the SG stars as the gaseous ejecta either 
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by massive asymptotic giant branch (AGB) stars ( "AGB 
scenario" ; iVentura et al.llMlt iD^Antona fc Ca"ioi|[200l 
iKarakas et alj 120061) or by fast rotating massive 
stars ("FRMS scenario" ; Ea ntzos & Charbonncl 2006:; 
iMevnet et alll2006t iDecressin et al.il2007i) . These stud- 
ies investigated whether the observed physical properties 
of the Galactic GCs are consistently accounted for in 
the models. Problems ar e present in either sc enario (for 
a critical discussion, see lD'Ercole et alll2010[) . although 
some chemical signatures of the SG, such as the Lithium 
abundances (jD'Orazi & Marino 20K); Shcn ct al. 2010) 
or the magnesium depletion (IVentura fc 'D'Antonall2009i : 



IVentura. Carini fc D'Antonali2011[ ) favour the AGB sce- 
nario, in which the hot CNO processing occurs by "hot 
bottom burning" (HBB) in the convective envelopes, 
that are recycled to the intra-cluster medium by wind 
mass loss. The AGB ejecta are likely to be converted 
into new stars, thanks both to their low ejection veloci- 
ties and to the fact that they evolve after the end of the 
energetic Supernovae II epoch. 

Most GCs sho w none or sca r ce iron spread among 
their stars (e.g. ICarretta et al.l l2009cl ). A few clus- 
ters, however, must be born in a more complex "pro- 
tocluster" environment. It has been often suggested 
that uj Gen, the most massive galactic GG, formed 
in a now dispersed dwarf galaxy (e.g. iNorris et al.l 
fl997t IRev et al.l [200l IBekki fc Noi^ I2006D and, pos- 
sibly, has an evolutionary history similar to that of 
M54, whose host galaxy, the Sagittarius dwarf, is now 
merging into the Galaxy. Both lo Gen and M 54 
show a large spread in iron, a; Gen contains stars 
with iron contents spanning a range from [Fe /H] 
2 to [Fe/H1^ 



-0.6 (Norris ct al. 1996; Suntzcff fc Kraft] 
1996'; 'Lcc ct aHlggQ; Pancino ct al. 2000; Sollima ct alj 
2005a: Johnson ct al. 2008; Johnson fc Pilachowski 201(| 
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The complexity of uj Gen star formation episodes 
is also shown f rom its multipl e red giant and sub 
giant bran ches f Lee et al.l 119991: lPa.ncino et al.l 120001: 
iSolhma et al. 2005bD. iBedin et al.l (|2004D found that 
the main sequence (MS) contains a well separated blue 
bra nch, probably then hiding a large overabundance of 
He (INorrisll2Q0l . as its stars are more iron-rich than 
stars on the redder side (jPiotto et al.ll2005l ). A third, less 
popul ated MS (MSa) has been discovered bv IBedin et al.l 
(|2004l ) on the red side of the red MS and it has bee n asso- 
ciate d to the most metal rich stars (e.g. Villanov a et al.l 
|2007| ). As well as blue MS stars, the colors and the metal- 
licity of the MSa suggest that it is populated by He rich 
stars ()Bellini et al.ll2010D . 

The presence of large variations in the light elements 
abundanc es and s-process elements is well documented 
see, e.g. iPaltoelou fc Norr is' '1989'; 'Norris fc Da GosT" * 



20001 ). Norris fc Da Costal (|199. 



[19951: ISmith et al 
show that w Gen giants also display huge aluminum 
variations, anticorrelated with oxygen and correlated 
with sodium, another typical pattern present in mono- 
metallic GCs. More rece ntly, the spe ctral analysis of a 
sample of 840 gia nts bv iJohnson fc Pilac howski (2010) 
and 300 giants by[M arino et al.l (|201l] ) shows with more 
clarity a very interesting peculiarity. Although the pres- 
ence of a wide range of iron contents points towards a 
chemical evolution similar to that of galactic environ- 
ments, the stars in each bin of iron content for [Fe/H]<- 



1.3 also show the typical 0-Na anticorrelation of GCs 
having a homogeneous iron content. On the contrary, at 
[Fe/H]>-1.3 the O and Na abundances show a direct cor- 
relation, and the [Na/Fe] values are always much larger 
than the solar ratio. 

The attempts to model the chemical evolution of 
UJ Gen, in the framework of a dwarf galaxy environment , 
have not been fu lly successful (e.g. iRomano et all 120071 : 
iChoi fc Yil [2008D especially for the difficulty in model- 
ing the origin of the very helium-rich population(s). The 
most recent modelling by Romano et al. (2010) docs not 
yet include the spectroscopic results concerning the O- 
Na anticorrelation for each subset of iron content we 
just discussed. In fact, a linear chemical evolution path 
for each element can not describe this complex situa- 
tion. Their work has to make a series of complex hy- 
potheses concerning winds that must preferentially eject 
some products of nuclear evolution, in order to reproduce 
some characteristics of the cluster stars. In particular, a 
timescale of a few Gyrs is needed to form the stellar gen- 
erations with different iron content in this cluster, and 
this rules out the massive AGB scenario for the formation 
of the super-helium rich popula tion(s); the formation 
of this population, according to IRomano et all ()2010l ) 
should t hen be a t tribut ed to the yields of FRMS, but 
see, e.g., iRenzinil (|2008[ ) against this possibility. Since 
our work is focused on the AGB scenario also for the 
formation of the very helium rich population(s), we will 
dis cuss later on ( SectjSt why this fails in the approach 
bv IRomano et al.l (j201Q) . 

We point out another ingredient of the w Gen puzzle: 
several studies have shown that the blue MS stars are 
more strongly concentrated in the cluster core than the 
red MS stars, and that this higher concentration holds for 
the most metal rich populations too (fSuntzc fF fc Kra ^ 
1996; Nor ris et al.l r i996: Rev ct al. 20041 iBellini et all 
2009; Johnson fc Pilachowski, 2010. '). These observed dif 



ferences in the spatial distribution of first- and second- 
generation stars are consiste nt with the p redictions of the 
simul ations presente d in D' Ercole et al.l ([2008). Specifi- 
cally, iD'Ercole et"all ([20081) show that the SG formation 
in "simple", mono-metallic GCs can occur in a cooling 
flow that gat hers the mas sive AGB ejecta in the cluster 
core (see also lBekkill2010l for 3D simulations confirming 
t he results of D'Ercole et al. 2008). 

ID'Ercole et al.l (|2010l ) extended the dynamical model 
to predict the abundances in the SG and, in particular, 
how the 0-Na and Mg-Al anticorrelations are built up. 

The remarkable 0-Na-Fe abundance distribution de- 
scribed above, introduces a new complex ingredient in 
the evolution of uj Gen and suggests that a new approach 
is needed to model the chemical evolution of this cluster. 
We suggest that, despite the complexity of the abun- 
dance patterns in uj Gen stars, the formation history 
of this clus ter may be exp l ained according to the key 
steps of the ID'Ercole et al.l ()2008[ ) model, that also pro- 
vides the key to understand how the stars with [Fe/H]>- 
1.3 may be directly born from the pure eje cta of mas- 
sive A GB stars (as already recognized by Johnson et al.l 
|2009| ). We also show that the cluster formation outlined 
in this work is in agreement with the abundances of other 
light elements (Mg and Al). It is also consistent with 
the spatial distribution of uj Gen populations. The criti- 
cal problem of the formation timescale(s) of the multiple 
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Fig. 1. — Literature data for O, Na, Mg and Al as a function of the iron content [Fe/H], The lower boundaries for Na and Al and the 
upper boundaries for O and Mg both follow closely to the trends predicted from SN II ejecta, weighted by a "typical" initial mass function. 
The top left panel shows the Mg data by Norris & Da Costa (1995), the bottom left panel shows the Na data by Johnson & Pilachowski 
(2010) (triangles) plus the data by Marino et al. (2011) (full squares). The top right panel shows the Al data by Johnson et al. (2008, 
2009) (cyan triangles) and Johnson & Pilachowski (2010) (green triangles), and the bottom right panel shows the Johnson (2010) O data 
(triangles) plus th e data by M arino et al. (2011) (full squares). The full big dots (red) in the panels represent the AGB yields as a function 
of [Fe/H] from I Ventura fc D'Antonal i2009i) . and in particular the minimum Mg and O yields, and the maximum Al yield. For sodium, we 
plot two values: the yields for SMq (lower line) and 4.5 M©. 

populations constituting the GC as we see it now is fi- 
nally discussed in SectlSj where we remark that the abun- 
dances of s-process elements as a function of metallicity 
may be uncorrectly biasing our views. 

2. THE P-CAPTURE DATA IN ui Cen 

While in this work we will focus mainly on the Na 
and O abundances in uj Cen, and will use the data 
from iMarino et al.l ()2011l ) for the sake of homogene- 
ity, here we provide some more information on re- 
sults for Mg and Al, whose abundances constitute an- 
other powerful indicator of p-capture processing in the 
gas from which the SG stars formed. Fig. [T] shows 
the O, Na, Mg and Al abundances as a function of 
[Fe/H] from different data, and the lines corresponding 
to theoretical models. These latter will be di scussed in 
Sect. [4.11 Mg data are tak en from iNorris fc Da Costa 
(119951). and Al d ata are from JJohnson et al.l (|2008ll2009l) : 
iJohnson fc Pilachowski (201Q| ). The figure shows that Na 
and O data from i Johnson fc Pilachowski (2010) are qual- 
itatively compatible with the .Marino et al- (,201 li) data, 
although a trend to have lower N a and O value s is ev - 
ident, especially at [Fe/H]>-1.3. IMarino et all (|201lD 
show in their Fig. 6 that t heir own oxygen dete r mina- 
tion s can be larger than the IJohnson et al.l (|2008l . I2009( ) 
and IJohnson fc Pilachowskil (|2010( ) determinations, and 
the difference may reach 0.3 dex or more at low oxygen 
values. Similar systematic (although smaller) differences 
appear also in the sodium determinations. For this rea- 
son we decided not to combine the different datasets, 
and we mainly use our own data to put constraints on 
the 0-Na (anti)correlation(s). 

We show the and Na abundances from the data 
by IMarino et al.l (|201lD in Fig. [5] We consider 6 bins 
of metallicity, keeping in mind that the errors on the 
abundance determination may easily shift some points 



from one bin to the adjacent one(s) ([Marino et al.l[20Tll ). 
We plot [0/H] and [Na/H] instead of the values nor- 
malized to the iron abundance, scaling all the values to 
the lowest metallicity [Fe/H] =-2. In order to increase 
the rel atively few data at [Fe / H]>- 1, we add the val- 
ues by IJohnson fc Pilachowskil ([2010f ) for these metallic- 
ities. These data increase the scatter, but, in this kind 
of plot, they do not change the qualitative behaviour of 
the metal rich sub-populations. We point out here some 
interesting patterns emerging from this figure (see also 
the right panel of Fig. |3|): 1) the lowest metallicity bin 
shows a scarce, if any, anticorrelation; 2) each interme- 
diate metallicity bin in the range -1.9 ^ [Fe/H]<-1.3 
shows independent extended 0-Na anticorrelations; 3) 
the highest metallicity bins (apart from a few stars in 
the range -1.3 ^ [Fe/H] < -1.1 located in the 0-rich, 
low-Na region) show very clearly a radically different be- 
haviour, a direct 0-Na correlation, with oxygen abun- 
dances slightly reduced and very large sodium values, up 
to ^20 times the solar ratio. 

The huge sodium values of the most iron rich bins is 
a very atypical feature; we are not aware of any other 
galactic environment in which such a feature is found. In 
particular, the field stars of the Sagittarius dwarf galaxy 
having comparable iron contents, have solar, or lower 
than s olar, [Na/Fe] ([Sbordone et al.ll2007l:lCarretta et al.l 
l2010a| ). Thus the field of the Sagittarius formed in a 
very different context than the metal-rich population 
in w Cen, although other elemental ab undances seem 
to sug gest a similar chemical evolution ([Carretta et alj 
I2010bf) . 

In Fig. [2| we also plot the yields of the massive AGB 
models (see Fig. |3|for the yields down to ^Mq), scaled 
in the same way as the data. We limit the masses to the 
range 4.5<M/ Mq<6.5. Both O and Na increase when 
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Fig. 2.— The Na vs. O data by Marino et al. (2011) for the 
giants in ui Cen are subdivided into different iron content groups, 
and plotted with the abundances are normalized to [Fe/H]=— 2. 
The open triangles are the data for [Fc/H]>-1 (large triangles 
from Marino et al. (2011,), small ones from Johnson & Pilachowskil 
pOlO)). The lines with symbols on the left arc the yields from the 
models by Ventura and D'Antona (2009) for the labelled metallic- 
ities. Masses from top to bottom of each segment are 4.5, 5, 5.5, 
6Mq. The last point is 6.4 Mq (Z=6x 10"''), 6.5 AfQ(Z=xlO-3) 
or 6.3M0(Z=4xlO-3). 

the stellar mass decreases (see Fig. [3] for the super- AGB 
behaviour). In fact, smaller masses have lower temper- 
ature at the HBB boundary of convection (Tbce), and 
only allow a milder ON and Na cycling. Increasing the 
metallicity, the models define parallel and partially su- 
perimposed diagonal lines, that shift to larger O and Na 
values for larger metallicity. This is due to two con- 
comitant reasons. The first reason is that, by increasing 
the metal content, the average content of each element 
increases. The sodium yield is a result of four differ- 
ent processes occurring in the HBB envelopes: 1) second 
dredge up of the sodium produced in the stellar inte- 
rior in the previous phases of evolution; 2) burning to 
sodium of the ^^Ne achieved during the same process of 
second dredge up; 3) and burning during the HBB phase 
along the Ne-Na cycle. 4) In addition, in the relatively 
small masses (M<4.5Mo), sodium is also produced by 
HBB of the ^^Ne dredged in the envelope at each third 
dredge up process. As sodium is produced from ^"^Ne 
and Ne, and this latter is ultimately produced from 
^''N, the average production of ^^Na, for similar HBB 
conditions must increase with stellar metallicity. For the 
same reason, the equilibrium value of O in the full CNO 
processing is larger at larger C, N and O abundances. 

The second reason is that the HBB efficiency decreases 
with increasing metallicity, as the higher opacities affect 
the envelope structure and allow smaller Tbce's, lower 
ON p rocessing and Na depletion () Ventura fc D'Antonal 
|2009| ) . Notice that the computation of the sodium yields 
is a delicate issue, as they depend on uncertain cross 
sections, and also on the initial abun dance of ^°Ne. 
This p roblem is examined extensively by iD'Ercole et al.l 
(|201Cll) . who adopt Na abundances larger by 0.2-0.3dex 
with respect to those plotted in Fig.|3l Yields from other 
modelers (e.g. iKarakas fc LattanziQ,2007, ) attain smaller 



oxygen depletion and larger sodium abundances for a 
fixed mass, due to the different assumptions made for the 
convective e fficiency and mass los s . This issue is largely 
discussed in lVentura fc D'Antonal (|2005al lbl). 

The theoretical yields seem to define a border limit to 
the abundances of the data, and the yields of the mas- 
sive AGBs of the highest available metallicity Z=4 x 10^^ 
go through the data in the two bins of largest iron con- 
tent. As the model yields reproduce the observed 0-Na 
correlation for the high-metallicity population, we are 
led to suggest that the stars at [Fe/H]>-1.3 form from 
the pure ejecta of massive AGBs of the same metallic- 
ity^ and this may represent a key ingredient to under- 
stand the key l ines of chemical evolution in w Cen. If 
we had adopted [Johnson fc Pilachowskil (|201ClD data, the 
same conclusion would have been reached, but only for 
[Fe/H]>-1, as their oxygen values for giants having - 
1.3<[Fe/H]<-l can reach much lower abundances, likely 
due to systematics both in [Fe/H] and fO/Fe] betwee n 
the two data-sets (see Figs. 4-6 in lMarino et al.ll201lD . 
We now try to use this finding as a constraint for the 
formation model of the cluster. 

3. THE O-NA PLANE: DATA POINTS, THEORETICAL 
YIELDS AND PREDICTIONS OF D'ERCOLE ET AL. 
MODEL 

Before we discuss a possible modelling of the 0-Na 
data in uj Cen, it is propaedeutic to summarize the status 
of art of modelling the 0-Na anticorrelation in single- 
m etallicity GCs. The left panel of Fig. |3| shows the data 
bv lCarretta et al.l (j2009aD for giants in these clusters. Su- 
perimposed, we plo t the whole oxygen and sodiu m yields 
of AGB stars from 'Ven tura fc D'Antonal (|2009 '). As dis- 
cussed in Sect. [21 masses increase from top right to bot- 
tom left. For the Z=10~^ m odels, we also plot the results 
of the yield computation bv lVentura fc D'Antonal (|201lD 
for the masses 6.5<M/ Mq<8, that, after having ignited 
Carbon in conditions of semi-degeneracy, and formed an 
0-Ne degenerate core, follow the super-AGB phase. We 
see that these more massive progenitors do not change 
the location of the Na-0 model correlation, but simply 
reverse it, and now the O and Na content increase for 
increasing mass. This is due to the huge mass loss rates 
of super-AGBs. So by increasing the initial mass from 
3 Mq to the maximum mass that does not explode as 
supernova, the Na-O line first goes down, reaches a min- 
imum value in the yields (at masses ~6.5 Mq), and then 
reverts back to larger O and Na. 

Three main features are evident from the compari- 
son between yields and single-metallicity globular cluster 
data: 

1. the AGB models of different masses provide a di- 
rect correlation between sodium and oxygen abun- 
dance in their ejecta. The presence of the O- 
Na anticorrelation in the GC data is generally at- 
tributed to dilution of the hot-CNO processed gas 
of the stellar ejecta with "pri stine" gas ( Bckki et al] 
200?; 'D' Antona fc Ventural 120071 : [P'Ercole et ah! 
2008. 20l'oi r Dilution is also needed when dealing 
with FRMS ejecta (iPrantzos fc Charbon^ 120061: 
iDecressin et~aLll2007l: iLind et al.ll201lD . In ah clus- 
ters dilution seems to be present. Detailed hydro- 
dynamical models able to explain the dynamics and 
the presence of pristine gas during the SG forma- 
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Fig. 3. — Left panel: Na vs. O data for GCs in Carretta et al. (2009a) are shown as open triangles. Superimposed, we plot the yields 
for AGBs from Ventura & D'Antona (2009). Metallicities are labelled. Masses from top to bottom of each segment are 3, 3.5, 4, 4.5, 5, 
5.5, 6AfQ. The last point is 6.4 Afg (Z=6x 10"*), 6.5 M0(Z= X lO'^) or 6.3 Mq (Z=4x 10-^). The Z=10-^ yields are completed with the 
models by Ventura & D'Antona (2011) for super— AGB models of 6.5, 7, 7.5 and SMq, that revert upwards to larger oxygen and sodium 
abundances. Right panel: The Na vs. O data by Marino et al. (2011) for the giants in lo Cen are plotted on the data of mono-metallic 
GCs shown in the left panel. The stars are subdivided into different iron content groups. The l i ne rep resents the chemical pattern of the 
SG in a monometallic cluster with Z=10~'^, according to the standard model bv .D'Ercole et al.l II2008I) . updated using the model yields of 
the left panel and the prescriptions described in the text. Dots label ages of 80, 100 and ISOMyr from the formation of the FG. The line 
ends at 200Myr. 



tion of all clusters ar e still lacking. This issue is 
discussed in detail bv lD'Ercole et al.l (|2011[ ). 

the yields for masses M ;S 4.5-5 Molie outside the 
plane region occupied by the 0-Na observational 
values. This represents a powerful constraint of the 
models, showing that the winds of massive AGB, 
for some reason, can not contribute to form SG 
stars at the epoch in which these -or smaller- 
masses begin to evolve (>130Myr). A powerful 
similar constraint is imposed also by the observed 
constancy of the sum of CNO abund ances among 
the stars of many GCs (e.g. M92 (iPilachowski^ 
1981 . NGC 288 and NGC 862 (IDickens et aL 
1991. M3 and M13 (|Smith et al.l I1996D . M4 
(jlvans iEaD[i999)) The occurrence of numerous 
third dredge up events during the AGB evolution 
of masses M<4.5 M p, (Ventura fc D'An tona 2009; 
iD'Ercole et"aIll2010D increases the total CNO. 

the yields do not represent well the lowest oxy- 
gen val ues th at can be present in the data. 
iD'Ercol e et al.l ([2010;) assumed semi empirical 
yields down to [0/Fe]= —1 in the super- AGB 
masses, by adopting a 'deep -mixing" suggestion. 
iD'Antona fc Ventural ()2007l ) attribute the very 
small oxygen abundance present in the giant stars 
of the most extreme SG to deep extra-mixing act- 
ing in the progenitors having a very high helium 
content. In these stars, the small molecular weight 
discontinuity during the red giant branch evolution 
may not be able to preclude deep mixing and low- 
ering of oxygen in the envelope. Stars having the 
largest helium abundances would then populate the 
blue MS and show the strongest chemical anoma- 



lies (see later). The presence of [0/Fe] abundances 
smaller than those predicted by AGB models only 
among the giants, and not amo ng the turnoff 
or su bgiant stars so far examined (jCarretta et al.l 
120061 ) , and in stars belonging to clusters also show- 
ing the signature of a very helium rich MS, such 
as NGC 2808 and uj Cen, are a possible argument 
in favour of this interpretation. Obviously, the de- 
termination of a very low oxygen content — similar 
to the one in the most anomalous red giants — in 
the atmosphere of blue MS stars in these clusters 
would falsify this suggestion ()Bragaglia et al.|[201Cll 
were not able to measure oxygen in their spectrum 
of the blue MS in NGC 2808.). 

Let us discuss how the 0-Na ant icorrclation is intro- 
I (|2008ll2010l ). In order to explain 



duced in lD'Ercole et al.l 
the data of NGC 2808, they propose a model in which 
the "extreme" SG, corresponding to the most oxygen- 
poor stars and to the stars populating the very helium 
rich blue MS, are born directly from the pure ejecta of 
the super-AGBs of the FG. In fact, after the SN II epoch, 
that has cleared the cluster from the residual g cool- 
ing flow sets in and brings to the core the low velocity 
winds of the super-AGBs, so that stars can form directly 
from their very helium rich and hot-CNO processed ma- 
terial. Then a phase of mixing of the massive AGB ejecta 
with pristine gas (reaccreted in the cluster core) follows, 
giving origin to the "intermediate" SG stars, with milder 
oxygen depletion. The SG formation stops due to the 
onset of the Type la SN epoch before the pristine gas 
is exhausted. In the right panel of Fig. |31 the slanted 
V-shaped line s hows an example of th e 0-Na patterns 
predicted by the ID'Ercole et al.l ()201Cll ) model, based on 
the whole set of yields for masses from 3 to 8 Mq shown 
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in the left panel, and extended for 200Myr after the FG 
formation (D'Ercole et al. 2011, in preparation). Time 
along the line proceeds from the left side to down right, 
and then climbs up towards the upper right side. The 
"extreme" SG is the left side of the V-shaped line, up 
to the point at t=80 Myr. The "intermediate" SG stars 
follow the line from 80 to 100 Myr. If star formation pro- 
ceeds further, the subsequent SG would follow the right 
side of the V, from t=100 Myr to the end, when stars 
form directly from the ejecta of the evolving AGBs, so 
that the simulation follows the direct correlation 0-Na 
shown in the left panel, but this side of the curve has no 
observational counterpart in "normal" GCs. The simu- 
lation stops at 2xl0^yr, when the evolving AGB mass is 
~3.7M0. The mo del is plotted on the [0/Fe], [Na/Fe] 
data for cj Cen from lMarino et all (j2011|) , and we see that 
data points are present also on the upper side of the sim- 
ulation curve, but they have larger metallicity than stars 
showing the anticorrelation(s). This suggestive compar- 
ison requires some deeper understanding. 

4. INTERPRETATION OF THE uj Cen DATA IN TERMS OF 
D'ERCOLE ET AL. MODEL 

We discussed in the previous Sect. |3] how the 0-Na 
anticorrelation in the GC data is attributed to dilution 
of the hot-CNO processed gas of the stellar ejecta with 
"pristine" gas. Where does this pristine gas come from? 
iD'Ercole et aD (|2008| ) suggest that it may survive m a 
torus that collapses back on the cluster after the SN II 
epoch. If the system is much more massive than a typ- 
ical proto-cluster as in the case of ui Cen (possibly in- 
cluding a dark matter halo) the SN II ejecta (and type 
la SN ejecta as well) can not be lost outside the poten- 
tial well o f the pro t oclus ter. The 3D hydro simulations 
by Marcol ini et al.l (|2006[ ) show in fact that the collapse 
back includes the matter enriched by the SN II ejecta. 
Their model was first developed to be compared with 
Draco's dwarf galaxy data. Afterwards, they attempted 
to m odel directly the case of uj Cen (|Marcolini et al.l 
l2007f l. Star formation in the cluster region stops when 
SN II explode, and a successive burst occurs when the 
interstellar medium (enriched by SN ejecta) flows back 
down the cluster potential well. Thus the gas content 
in the cluster center has an oscillating temporal pro- 
file whose period is given by the time-interval between 
two successive starbursts. Each star burst provides star 
formation with increasing metallicity. If we add to this 
model the possible role of the winds of the AGBs evolv- 
ing in the progenitor systems, these may account for star 
formation of the very helium rich stars when the clus- 
ter core is still devoid of matter (before each collapse 
back of the ISM), and for the formation of stars from 
AGB matter diluted with the ISM enriched by the SN II 
eje cta, similar to the eve nts in the one-burst only model 
bv lD'Ercole etall (|2Q10f ). Suppose now that, after a last 
burst of star formation, there is no longer gas that can 
collapse back from the outskirts of the cluster. Never- 
theless, there can be further episodes of star formation 
thanks to the cooling flows of the AGB winds. Whatever 
is the detailed mechanism by which the 0-Na anticorre- 
lation(s) are built in uj Cen for [Fe/H]<-1.3, the 0-Na 
correlation we see at larger metallicities is telling us that 
the ISM in the outskirt of the cluster must be so reduced 
that only the AGB winds take part in the cooling flow 



and that these winds belong to the massive AGBs of the 
more metal rich populations. This interpretation only 
requires that these AGBs are delayed until there is no 
longer any pristine gas available to collapse back in the 
core of u! Cen. A possible hypothesis is that the grav- 
itational iiitera£twn_oftl^ Cen progenitor with the 
Galaxy (jDinescu et al.llT999l ) caused the loss of a signif- 
icant fraction of its mass, or of the dark matter halo, so 
that the last SN II burst, or the type la SN explosions 
could clean completely the central region and the whole 
proto-cluster from its residual gas. 

Following this scenario, we now can re-interpret the 
path of the simulation shown in Fig. |3]and explain why it 
resembles the shape of the u! Cen data in the Na-0 plane. 
In a mono-metallic cluster, the stars on the upper branch 
would be those formed by smaller mass AGB stars ejecta 
(M ^ 4.5 Mq), evolving at a later time, when the diluting 
gas is over. In ui Cen, the high Na abundances belong to 
the stars in the most metal rich bins, and correspond to 
the yields of (again) more massive AGB stars for those 
metallicities (Fig. [2]). The only requirements to explain 
why the most iron-rich giants in uj Cen have large sodium 
and moderately low oxygen are: 1) these stars are formed 
in a cooling flow directly from the ejecta of the massive 
AGBs having the same metallicities; 2) the AGBs that 
provide the ejecta evolve at the time when all the gas 
surrounding the cluster has been already exhausted. 

4.1. Consistency of the scenario: the magnesium and 
aluminum abundances 

We have so far focused our attention on the Na-0 
(anti)correlation(s) to outline a possible scenario for the 
formation of the multiple populations in uj Cen because 
of the wealth of data available for the abundances of these 
two elements. 

In this section, we show that existing data on other 
light elements lend further support to the proposed sce- 
nario. In Fig. [1] the top panels show the Mg and 
Al data as a function of [Fe/H]. We plot on the data 
the mini mum Mg yield and ma ximum Al yields from 
the Ventura fc D'Antonal (|2009[ ) models. These val- 
ues then represent the maximum degree of Mg burning 
and Al production i n massive AGBs of each metallic- 
ity. The Mg data by iNorris fc Da Costal (|1995[ ) show a 
remarkable constancy of magnesium for all metal con- 
tents, but the few low Mg points are at iron con- 
tent for which models predict Mg depletion. Quanti- 
tat ively, the models are a facto r ■-^2 o ff the data, but 
see lVentura. Carini fc D'Antonal ()2011h for a discussion 
of the dependence of the Mg depletion in AGB and 
super-AGBs on the input physics of the models. It 
is important to point out that no Mg depletion is pre- 
dicted at the highest metallicities, for which, according 
to our model, the SG formation occurred directly from 
the pure AGB ejecta. The upper boundary of Al data is 
in good agreement with the model yields. In particular, 
a smaller Al enhancement is predicted at high metallic- 
ity than at the lower ones, as observed. The Al data are 
also consistent with the hypothesis that some dilution of 
AGB ejecta with pristine gas occurs up to [Fe/H] ;S -1.3, 
while either a dichotomy in the Al values is found at - 
1.3 [Fe/H] ^ -1, or values compatible with pure ejecta 
at [Fe/H]>-1. 
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5. THE TIMESCALE OF FORMATION OF THE u) Cen 
POPULATIONS 

We do not attach a fixed timescale to the model pre- 
sented to account for the high sodium of the metal rich 
populations in lo Cen. Its main requirement is that the 
time delay for the formation of the massive AGB progen- 
itors in these iron bins must be of the order of the time 
required for the preceding formation events to exhaust 
the galactic gas in the neighborhood of the cluster, so 
that this gas does not collapse back to the cluster core 
after the last SN II epoch. We think that, in order to 
satisfy these requirements, the timescale for the forma- 
tion of all the populations in lo Cen must not exceed a 
few lO^yr. 

Notice that the upper boundary for [0/Fe] (see the 
right panel of Fig. [3]) increases with increasing metal- 
licity, in spite of the presence of the many 0-poor gi- 
ants defining the 0-Na anticorrelation(s). Standard 
galactic chemical evolution shows that [0/Fe] remains 
> until the iron is produced by SN II, and shifts 
to the solar ratio only when Type la SN, produc- 
ing much more i ron and no oxygen, beco me the main 
source of metals (jMatteucci fc GreggiolfT986ll . Thus the 
iron chemical evolution has been dominated by Type 
II SN, suggesting that the initial chemical evolution in- 
side the progenitor system occurred on a fast timescale. 
As also the other a-elements seem to go towards the 
solar ratio only for the most metal ri ch population, 
corres ponding to t he MSa and RGBa ijPancino et al.l 
120021: lOriglia et al.l [20031 ) . the timescale for the forma- 
tion of the last, metal richest, population in lo Cen 
must be in fact very close to the timescale for the on- 
set of SN la explosions. Although this timescale is 
uncertain, its value is obser vationally estimated to be 
in the range 40-300 Mvr (iMatteucci fc Recchil 120011: 
iMannucci et al.l 120061: iRaskin et al.ll2009D. and certainly 
not longer than SO OMyr (e.g. IGal- Yam fc Maod 12004 
ISullivan et al.|[2006[ ). Simulations adopting the double- 
degenerate model for this event show that the peak of 
the d istribution is reached after only 250Myr (Tornambe. 
|1989() . It is difficult to see how this epoch can be de- 
layed by more than a few 10* yr in lo Cen. In the 
iMarcolini et al.l (|2006f ) and in other similar models, the 
presence of a dark matter halo prevents the clearing of 
gas by type la SN, so some other process must have pre- 
cluded further star formation in the proto-w Cen (and in 
spheroidal dwarf galaxies as well), e.g. the gravitational 
interaction with the Galaxy. 

As the iron production is delayed in time and massive 
AGBs of different metallicity evolve at different epochs, 
also less massive (M<4.5M0) AGBs of the lower metal- 
licity bins would contribute to the cooling flow. These 
may be, e.g., some of the stars in the bin -1.3<[Fe/H]<- 
1). To understand the role of smaller AGB masses of 
lower metallicities, notice that the metal content of the 
massive AGB of high metallicity dominates over the 
anomalies of these latter stars (see the discussion of the 
0-Na yields shown in Fig. [2]). As a net result, the [0/Fe] 
and [Na/Fe] of the gas may be not affected too much by 
the yields of the less massive stars and the latter contri- 
bution to the average abundances remains a small per- 
turbation. Of course, a full model is needed to quantify 
the extent of this problem. 



HR diagram d ating resu lts in age differences from less 
than 2 Gyr fe. g.lLee et al.|[20 05: SoU ima et al.ll2005al) to 
5 Gyr or more (jVillanova et al.ll2007i) in recent literature. 
The influence of different helium content and/or of the 
total C-|-N-f O content on the isochrones turnoff level may 
easily change these age determinations (iD'Antona et al.l 
l200a iPietrinferni et all [20091) and bring them to much 
closer values. Joo, S. & Lee, Y.-W. (2011, in prepara- 
tion) have indeed shown that it is possible that all the 
populations are coeval, by taking into account the com- 
bined effect of helium and metallicity. The result is valid 
within the errors of ±0.3 Gyr that affect the determi- 
nation of relative ages from the turnofFs' location. This 
require s that the MS - a is in deed helium rich, as suggested 
also bv lBellini et"all (|20lol ). 

We address now the most important clue, that 
we consider a red herring, for the timescale of for- 
mation of this peculiar cluster, namely the evidence 
of a very fast and large growth, by more than a 
decade, in the abun dances o f Lanthanum and B arium 
([Johnson k, Pilacho wski 2010 HMarino et al.|[20lll) in the 
very small range of metallicity between [Fe/II]~-1.9 and 
[Fe/H]'^-1.5. This very dramatic s-process increase with 
iron, at such a low metallicity, is not fo und in other galac- 
tic environments (iTolstov et al.ll2009D. a s shown in Fig- 
ure 17 of [Johnson fc Pilachowski! ( 2010[ ). [Marino et al.l 
([201 If ) suggest a rapid production of these elements. 
The r-process is excluded, as it would also affect Eu- 
ropium as well, while the [Eu/Fc] ratio do es not in- 
crease in the same meta l licity interval (Smith et al.ll2000l : 
[Johnson k, Pilachowskil [20100 . Another plausible pro- 
duction site are the cores of massive stars (weak s- 
process, e.g. [Raiteri et al.|[l993l and references therein), 
but this chain should mainly increase the lighter nu- 
clei in the s-chai n, like '^^Fe, ^'^Cu, or ^^Cu (see, e.g. 
iPumo et al.[[2010l ). while Cu in the metal rich stars of 
LO Cen remains at low abundances ([Cunha et al.[ [2002[ : 
[Pancino et al.l [20021 ). Further non-standard models for 
the production of these s-process elements may occur in 
rotating massive stars ^Pignatari et al. 2008). 

Nevertheless, the best acknowledged production site of 
these s-pro cess elements are t he lower mass AGB stars 
(M<3Mq) ([Busso et al.|[1999l ) that evolve on timescales 
of several hundreds of million years or on gigayears 
timescales, and these stars have been always considered 
the obvious responsible for this s-process production. Of 
course, if the -say- 3 Mq are responsible for the s-process 
production, this happens at a time at which all the pos- 
sible producers of helium and of the 0-Na anticorrela- 
tions (the stars of M^SM©!) have already evolved. 
We remark again that the 3 Mq stars would produce in- 
gent amounts of sodium, but do not deplete oxygen nor 
produce helium. It may well be that these masses con- 
tribute to increase the s-process abundances when the 
[s/Fe] value flattens (at [Fe/H] > - 1.5), but it is diffi- 
cult to envision how they can contribute to the extreme 
s-process increase that occurs at -1.5>[Fe/H]>-1.9, a 
range of metallicities that show independent, and even 
extreme, 0-Na anticorrelations (see the blue squares and 
the red triangles of Fig. [2|) and where there is no correla- 
tion of the s-process abundances with, e.g., the oxygen 
depletion. So the problem is that the timescale to build 
up the 0-Na anticorrelation(s) is much shorter than the 
timescale needed for standard s-process production sites. 
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We do not see an easy wa y out of this problem , and in 
fact this is the reason why iRomano et aLl ()201Cl[ ) had to 
exclude the AGB contribution for the formation of the 
helium rich population, and resort to the contribution of 
massive stars. 

We suggest that, in the same way as the 0-Na an- 
ticorrelation occurs only in GCs, and is due to the pe- 
culiar chemical evolution provoked by the collection of 
winds from a restricted range of initial stellar masses, 
we may be facing a process that has not yet been ex- 
plored in stellar models. A not well discussed site of 
s-nucleosynthesis occurs in the carbon burni ng shells of 
the t ail of lower mass progenitors of SNII fe.g. lThe et al.l 
l2007f) , and their contribution may peculiarly become ap- 
parent mainly in the evolution of the progenitor sys- 
tems of clusters like w Cen, that are unable to get rid 
of the SN II ejecta, but still give birth to a cluster in 
which cooling flows play a dominant role. We remark 
that the only similar trend of increase in the s-process 
abund ance, in the sa me ran ge of metallicities, is foun d in 
M 22 (: Marino et ah 2009; D a Costa fc: Marino|[2010l ) . al- 
though this cluster completely lacks the most metal rich 
populations ([Fe/H]>-1.5) present in lo Cen. Even the 
smaller cluste r NGC 1851 shows an s- process increase by 
a factor 2-3 (jYong fc Grundahill2008D . and this may be 
correlated to a po ssible, still discussed, s mall increase in 
the i ron content (jCarretta et al.|[2"010d : iVillanova et al.l 
120101) . It would be interesting to know whether the clus- 
ter NGC 2419, that harbors two popu lations with a small 
variation in Calcium content (.Cohen et al.ll2 0l"o|), and a 
popu lation with very high helium ((Di Crisc ienzo et al.l 
1201 It) also shows a similar s-process dichotomy. It looks 
very difficult that these very different clusters managed 
to evolve in such a similar way, if there is not an unify- 
ing model for the s-process production, linked mainly to 
stars that evolve into SN II. Further work is required to 
clarify this issue. 

6. CONCLUSIONS 

We have exam i ned th e patterns of the 0-Na-Fe data 
by iMarino et "aD (|201lD for the giants in w Cen. Most 
stars at the largest iron contents ([Fe/II]>-1.3) are all 
very sodium-rich and moderately oxygen depleted, a pe- 
culiarity not predicted in standard galactic evolution, 
and not found in other environments apparently simi- 
lar, like in the metal rich field stars of the Sagittarius 
dwarf galaxy. Furthermore, the metal-ric h stars pop- 
ulatin g the MS-a might be helium-rich (| Bellini et al.l 
I2010D . We have shown that these peculiar abundances 
(including high helium content) match those of the ejecta 
of the massive AG B stars of comparable metallicities 
(| Ventura fc D'Anto na 2009). 

We attribute the presence of a direct 0-Na correlation 
for the stars at [Fe/H]>-1.3 to formation, in a cooling 
flow, from the ejecta of massive AGB stars of the same 
metallicity, that evolve in the cluster when the gas in the 
surroundings has been all exhausted by previous episodes 
of star formation. These previous events have formed the 
populations with [Fe/II]<-1.3, from the ejecta of AGBs, 
diluted with the recollapsed gas, metal enriched by the 
SN II preceding burst. This process forms the typical 
Na-0 anticorrelation patterns. The reason why massive, 
metal rich AGBs are available for this phase of star for- 
mation is that, with respect to the AGBs at smaller iron 



content, they evolve with the time delay that has been 
necessary to form the metal-rich components in the pro- 
genitor system. Consequently, uj Cen differs from mono- 
metallic clusters because its initial mass must have been 
even larger (possibly including a dark matter halo) than 
that of a typical cluster. 

Consequently, uj Cen differs from mono-metallic clus- 
ters because its initial mass must have been even larger 
(possibly including a dark matter halo) than the ini- 
tial mass of a typical cluster. This led to an evolu- 
tion in which SN II ejecta could not be expelled from 
the cluster, but were incorporate d in successive bursts o f 
star formation, as described by iMarcolini et al.l ()2006D . 
The existence of helium-rich populations and the con- 
centrated spatial distribution of the helium-rich and of 
the most metal-rich populations (S untzcff & Kraft 199^ 
Norris et al.' 1996': 'Rev ct alj [20041: IBellini et al.l I200a 
Johnson fc Pilachow ski 201C|). however, shows that also 
in this cluster the gas from which second generation stars 
form collects in a cooling flow in the cluster central re- 
gions as described in the simulations by D'Ercole et al. 
(2008) for standard GCs. A further reason, to think 
that these latter stars were formed in cooling flows, is 
the 0-Na composition of their stars, that follows the di- 
rect 0-Na correlation predicted by massive AGB stellar 
models. We showed that the formation of the lo Cen 
mu ltiple populatio ns seems to follow the basic steps of 
the lD'Ercole et al.l (2010) model for mono-metallic GCs. 
The 0-Na patterns, and the Mg and Al abundances in 
UJ Cen then contribute to validate the AGB pollution 
model for the formation of multiple populations in GCs. 

An important and still unsolved issue in the forma- 
tion of ui Cen multiple populations is which are the stars 
responsible for the fast increase in the s-process abun- 
dances, that occurs at the epoch in which the metallicity 
raises from ~-1.9 to ~-1.5. We discussed the common 
view of s-elements production in low mass AGB stars, 
and concluded that it is in contrast with the other strong 
indicators of formation of the GC populations in AGB- 
dominated cooling flows. We suggest that the s-process 
patterns are a red herring, and propose that their evolu- 
tion is due to a process not yet identified or well explored 
in stellar models, such as, e.g., the nucleosynthesis in the 
carbon burning shells of the tail of lower mass progenitors 
of SNII. This (postulated) nucleosynthesis site becomes 
dominant in the chemistry of clusters showing a spread in 
metallicity, in the same way as the 0-Na anticorrelation, 
due to the peculiar chemical evolution provoked by the 
collection of winds from a restricted range of initial stel- 
lar masses, dominates the chemistry of GC stars. Further 
work is required to clarify this suggestion. For the same 
reasons, we suggest that the whole time of formation of 

Cen is at most a few lO^yr. 
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